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I. INTRODUCTION 

The g-state Potts model has served as a valuable model for the study of phase transitions 
and critical phenomena |jl|J^. On a lattice, or, more generally, on a (connected) graph G, at 
temperature T, this model is defined by the partition function 



ZiG,q,v) = J2e'^'' (1.1) 



with the (zero-field) Hamiltonian 



n = -jJ2s.^., (1.2) 



where Uj = 1, ..., g are the spin variables on each vertex i E G; P = {kBT)~^; and (ij) denotes 
pairs of adjacent vertices. The graph G = G{V, E) is defined by its vertex set V and its edge 
(bond) set E] we denote the number of vertices oi G as n = n{G) = \V\ and the number of 
edges of G as e{G) = \E\. We use the notation 

K = pj , a = e^=u-\ v = a-l (1.3) 

so that the physical ranges are (i) a > 1, i.e., v > corresponding to oo > T > for the 
Potts ferromagnet (FM), with J > 0, and (ii) < a < 1, i.e., — 1 < ^ < 0, corresponding to 
< T < oo for the Potts antiferromagnet (AFM), with J < 0. One defines the (reduced) 
free energy per site / = —(3E, where F is the actual free energy, via 

f{{G},q,v)= limln[Z(G,g,t;)V-] (1.4) 

where we use the symbol {G} to denote lim^^oo G for a given family of graphs. 

On a two-dimensional lattice, for the q = 2 Ising case, and for q = 3,4, the Potts ferromag- 
net exhibits a second-order phase transition from a paramagnetic (PM) high-temperature 
phase to a low-temperature phase with spontaneously broken symmetry and long-range fer- 
romagnetic order (magnetization). For g > 4, this transition is first-order, with a latent heat 
that increases monotonically with q, approaching a limiting constant as g — > oo 0,^. The 
behavior of the Potts antiferromagnet depends on the value of q and the type of lattice, as 
will be discussed further below. 

Let G' = {V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex 
set V and a subset of the edge set, E' C E. Then Z{G, q, v) can be written as the sum ||^ 

ZiG,q,v)= ^g'=(^\^(^') (1.5) 

G'CG 

where k{G') denotes the number of connected components of G'. Since we only consider 



connected graphs G, we have k{G) = 1. The formula (|1.5| ) enables one to generalize q from 
Z_|_ to R_|_ (keeping v in its physical range). The formula ( |1.5| ) shows that Z{G,q,v) is a 
polynomial in q and v (equivalently, a). The Potts model partition function on a graph G is 



1 



essentially equivalent to the Tutte polynomial 0- 0] and Whitney rank polynomial [Q], p- 
n| for this graph, and this connection will be useful below. 

Using the formula ( |1.5|) for Z{G,q,v), one can generalize q from Z+ not just to M+ but 
to C and v from its physical ferromagnetic and antiferromagnetic ranges < f < oo and 
— l<f<OtofGC. A subset of the zeros of Z in the two-complex dimensional space C^ 
defined by the pair of variables (g, v) form an accumulation set in the n — >■ oo limit, denoted 
B, which is the continuous locus of points where the free energy is nonanalytic. The program 
of studying statistical mechanical models with external field generalized from M to C was 
pioneered by Yang and Lee [O], and the corresponding generalization of the temperature 



from physical to complex values was initiated by Fisher [0. Here we allow both q and the 
temperature-like variable v to be complex. For a given value of v, one can consider this locus 
in the q plane, and we shall sometimes denote it as Bq, and similarly, for a given value of q 
(not necessarily G Z+), one can consider this locus in the plane of a complex-temperature 
variable such as v or 

It will be convenient to introduce polar coordinates, letting ( = |C|e*^. 

In this paper we shall present exact calculations of the Potts model partition function 
Z[G,q,v) for arbitrary q and v on self-dual strip graphs G of the square lattice with fixed 
width Ly and arbitrarily great length L^ with two types of boundary conditions. Letting 
L^ —>■ oo, we compute the resultant free energy and complex-temperature phase diagram. 
Results are analyzed for widths Ly = 1, 2, 3. We shall use these results to study the approach 
to the large-g limit of B. 

There are several motivations for this study. Clearly, new exact calculations of Potts 
model partition functions on lattice strips with arbitrarily large numbers of vertices are of 
value in their own right. This is especially the case since the free energy of the Potts model 
has never been calculated exactly for d > 2 except in the q = 2 Ising case in 2D. Just 
as the study of functions of a complex variable can yield a deeper understanding of func- 
tions of a real variable, so also the investigation of complex-temperature phase diagrams 
of spin models can provide further understanding of the physical behavior of these models. 
Besides [0, complex-temperature singularities were noticed in early series analyses (e.g., 
P^), and many studies have been carried out on complex-temperature (Fisher) zeros of 
the partition function of the Ising model and its generalization to the g-state Potts model 



1^ - [pl[| . In particular, several exact determinations of complex-temperature phase dia- 



grams of the Potts model on infinite-length, finite- width lattice strips pT| , [52| -p^,pq-|58| , in 
comparison with both exact solutions for the q = 2 2D complex-temperature phase dia- 
grams p!^ , |33|j37[] and finite-lattice calculations of Fisher zeros [^,^,^,^,0,^] have shown 



that, although the physical thermodynamic properties of these strips are essentially one- 
dimensional, one can nevertheless gain important insights into certain complex-temperature 
properties of the model on the corresponding two-dimensional lattice. For a model above its 
lower critical dimensionality, a complex-temperature phase diagram includes the complex- 
temperature extensions of the paramagnetic and ferromagnetic phases as well as a possible 



antiferromagnetic phase and other phases (denoted O in [33]) that have no overlap with any 
physical phase. For the infinite- length finite- width strips under consideration here, for finite 
g, the complex-temperature phase diagram includes only PM and O phases since there are 
no broken-symmetry phases. 

An early study of the complex-temperature phase diagram for the square-lattice Potts 
model led to the suggestion that the locus B, which is comprised of the circles |a ± 1| = \/2 



for g = 2 |T^ might generalize to the union of the circles | a — 1 1 = y/q and | a + 1 1 = -\/4 — q for 
1 < g < 4 p2|, but subsequent studies found that many of the Fisher zeros in the Re{a) < 
half-plane do not lie on a circular arc but instead show considerable scatter P5|P^P5|^ . 
The infinite square lattice is self-dual, and for calculations on finite sections of the square 
lattice, it was found to be useful to employ boundary conditions that preserve this self- 
duality. Stated more abstractly, one studies lattice graphs G with the property that the 
planar dual, G*, is isomorphic to G, which we write simply as G = G*. Here we recall that 
the dual of a planar graph G with n vertices, e edges, and / faces is the graph G* obtained 
by associating a vertex of G* with each face of G and connecting each pair of vertices on G* 
by edges running through the edges of G. It follows that n{G*) = f{G), e{G*) = e(G), and 
f{G*) = n{G). The Potts model partition function satisfies the relation 

Z(G, g, v) = gi-^(^)t;^(^)Z(G*, g, v,) (1.7) 

where the dual image of f is f^, given by 

Vd=- (1.8) 

V 

i.e., in terms of C, the duality map is the inversion map: 

= ^ ■ (1-9) 

Thus, it is also useful to plot Fisher zeros in terms of the variable C, since the accumulation 
set B is invariant under inversion for a self-dual graph: 

G = G* =^ B is invariant under C ^ - . (1-10) 

It was found that complex-temperature zeros calculated for finite sections of the square with 
duality-preserving boundary conditions (DBC's) have the appealing property of lying exactly 
on an arc of the unit circle |C| = 1 in the C, plane for Re{() ^ p9| , |38|j39| . (For physical 



temperature the coefficients of powers of a = e^ are positive, so there are no zeros on the 
positive real axis Re{a) > for any finite lattice. However, in the thermodynamic limit, the 
phase boundary crosses this axis at Oc = 1 + y^, i.e. Cc = 1-) In [^ several types of different 
self-dual boundary conditions were used in order to ascertain which features are common to 
each of these and hence might be relevant for the thermodynamic limit. However, because of 
the scatter of zeros in the Re{a) < half-plane and the dependence of the pattern of zeros, 
it has not so far been possible to reach a conclusion concerning the portion of the complex- 
temperature phase boundary B in this region (although some points on the boundary have 
been reliably located by analysis of series expansions ET 



An interesting exact result was obtained by Wu and collaborators, who gave an elegant 
proof [^, using Euler's identity for partitions, that for the Potts model on the square lattice. 



after having taken the thermodynamic limit so that zeros of Z in the complex C, plane have 
merged to form the continuous locus S, if one takes the further limit g — i> oo, this locus 
B is the circle |C| = 1 for the square lattice. It is straightforward to show that the same 
conclusion holds if one keeps Ly fixed and finite, and takes L^ — > oo, after which one takes the 
limit q -^ oo. Thus, for the infinite-length limit of the self-dual strips that we consider here, 
linig^oo B> is again the unit circle \C\ = 1. In 2D, the interior and exterior of this circle \(\ = 1 
form the complex-temperature extensions of the PM and FM phases, respectively. Since the 
infinite-length, finite-width strips under study here are quasi-one-dimensional systems, there 
cannot be any broken-symmetry phase at finite temperature for a spin model with short- 
range interactions and hence there is no FM phase or its complex-temperature extension. 
Note that any finite temperature point, i.e. — 1 < f < oo gets mapped to C = in the 
limit g ^ oo. The identification of phases thus proceeds as follows in this limit, for the 
infinite-length finite width strips: the interior of the circle \(\ = 1 is the PM phase, since 
it is analytically connected to the infinite-temperature point v = ( = 0. The phase in the 
exterior may be interpreted as equivalent to the T = point, in the sense that a finite value 
of ( in this region is obtained by taking the double limit q —>■ oo and |i;| ^ oo with v/^/q 
held fixed. 

Since, as noted above, for finite q, the actual pattern of zeros calculated on finite sections 



of the square in the Re{C) < half-plane show considerable scatter p9| , |38| , ^ , two questions 
arise naturally; first, having taken the 2D thermodynamic limit, if one starts with 1/q = 
and increases this quantity from zero, is there a finite interval in which the locus B continues 
to be the circle |C| = 1 before there are deviations, or do these deviations occur for any finite 
value of 1/q, no matter how small. We shall address this question here. A different question 
can also be posed for a finite section of the square lattice: for such a section, with a given 
size and given duality-preserving boundary conditions, is there a range in z = 1/q above zero 
in which all of the finitely many Fisher zeros still occur on the circle \(\ = 1 or not. This 
has been considered in PB| , P5| , |5(]| and we shall not pursue it here. Since one does not have 



an exact solution for the 2D Potts model for arbitrary q, and hence also no solution for the 
complex-temperature phase boundary B, it has not been possible to determine the precise 
behavior of this locus analytically in the q —>■ oo limit. 

Here one sees the value of exact solutions for the Potts model free energy and resultant 
complex-temperature boundary B on infinite-length, finite- width strips, since for these strips, 
one can obtain exact analytic answers to the behavior of B in the q ^ oo limit. As we 
shall discuss, it is easy to see that one aspect of this behavior is special to the quasi-one- 
dimensional nature of the infinite-length strips and is not relevant for the 2D model, namely 
that as 1/q increases from 0, a gap opens in the circle C = 1- This simply reflects the fact 
that for a quasi-one-dimensional spin model with short-range interactions there is no finite- 
temperature phase transition and the free energy is analytic for all finite temperatures, and 
hence for < C < cxd. However, this is not a drawback of the method, since finite- lattice 
calculations of Fisher zeros on sections of the square lattice have shown that they lie nicely 
on the circle |C| = 1 in the region near the point C = 1 (while avoiding the precise point C = 1 



if n is finite). Hence, one can infer tliat in the thermodynamic hmit, as 1/g is increased from 

0, this portion of B will remain as an arc of the unit circle. Indeed, from general arguments 

33| , one knows that for the model above its lower critical dimensionality, where there is 



a ferromagnetic phase, the portion of the phase boundary B that separates the complex- 
temperature extension of the paramagnetic phase from the FM phase must remain intact 
for finite as well as infinite q. If it were to bifurcate, this would imply a new third phase 
between the PM and FM phases, contrary to the known properties of the Potts model, so, 
given the invariance of B under the inversion symmetry ( p..lO|) , this portion must remain on 
the circle |C| = 1. 

Since the explicit calculations of Fisher zeros showed large scatter away from \C,\ = 1 in the 



Re{() < half-plane even for values well above g = 4, such as g = 10 ||39|j43[] , this suggested 
that the totality of Fisher zeros would only cluster on the circle |C| = 1 in the limit g — > oo 
itself but some zeros, and some portion of their accumulation set B would deviate from it 
for all finite q. Calculations using the usual Hamiltonian formulation of the Potts model and 
associated transfer matrix methods become increasingly cumbersome for large q because of 
the increasingly many states. For the purpose of studying the large-g behavior of the zeros, 
it is convenient to solve for Z{G,q,v) for arbitrary q and v on large finite lattice sections. 
This was done in [^] and the zeros were calculated for g up to 100; again these showed only 
a slow approach to the circle |C| = 1- Together with the previous calculations for q up to 10, 



it was concluded in HSlj] that this evidence supported the inference that the totality of Fisher 
zeros only lie on the circle \(\ = 1 in the limit g — *> oo. This type of calculations has also 
been done in [^ with the same conclusion. Our exact results on infinite-length finite-width 
strips complement these finite-lattice calculations and allow a rigorous conclusion that for 
these strips B deviates from \(\ = 1 for any 1/g no matter how small. 

These results are relevant in another way. Large-g series expansions (in the variable 
l/^/q) have been useful in studying the thermodynamic properties of Potts models ||62|- [^ . 



Large-g expansions (in the variable l/(g — 1)) have also been useful for a particular special 
case of the Potts model, namely the T = special case of the Potts antiferromagnet, where 
the partition function on a graph G reduces to the chromatic polynomial, 

Z{G,q,-l) = P{G,q) (1.11) 

where P{G, q) is the chromatic polynomial (in g) expressing the number of ways of col- 
oring the vertices of the graph G with g colors such that no two adjacent vertices have 
the same color [§,Q. Indeed, for a given graph G and for sufficiently large g, the Potts 
antiferromagnet exhibits nonzero ground state entropy (without frustration). This is equiv- 
alent to a ground state degeneracy per site (vertex), W > 1, since 5*0 = fc^lniy, where 
W{{G},q) = hm„^oo-P(G', g)^/". Large-g series expansions for Wr{{G}) = q^^W{{G},q) 
have been given, e.g., in [^- |]72| . 

A different motivation for the present study is the following. Just as was true for the 



double-complexification of field and temperature studied in [^- |^, where one gained a 
deeper understanding of the singular locus (continuous accumulation set of partition function 
zeros) by considering the separate projections in the planes of complex-field and complex 
temperature by considering these as parts of a single underlying submanifold (with possible 



singularities) in the C^ space of complex temperature and field, so also here, one gains similar 
insight into the projections of B in the complex q and v plane by relating these as different 
slices of the locus B in the C^ space defined by (g, v). 

II. GENERALITIES 

We refer the reader to our earlier papers containing exact calculations of Z{G,q,v) for 
a number of further details. We recall that the formal definition of the free energy may be 
insufficient to define this function at certain special special values q = Qs |^^ ; it is necessary 



to specify the order of the limits that one uses. We denote the limits with the two different 
orders as definitions using different orders of limits as /g„ and fnq'- 

fng{{G}, g, v) = lim lim n'^ In Z{G, g, v) (2.1) 

fqn{{G}, g, v) = lim lim n'^ In Z{G, g, v) . (2.2) 

Where necessary, we shall comment on which order of limits is used below. Of course in 
discussions of the usual g-state Potts model (with positive integer g) , one automatically uses 
the definition in eq. ( |1.1| ) with ( |1.2| ) and no issue of orders of limits arises, as it does in 
the Potts model with real g. As a consequence of the above noncommutativity, it follows 
that for the special set of points q = qs one must distinguish between (i) {Ba{{G}, qs))nq, the 
continuous accumulation set of the zeros of Z{G,q,v) obtained by first setting q = qs and 
then taking n ^ cxd, and (ii) {Ba{{G}, qs))qn, the continuous accumulation set of the zeros of 
Z{G, g, v) obtained by first taking n ^ oo, and then taking g — > g^. For these special points, 

Bnq ^ Bgn ■ (2.3) 

We have discussed this type of noncommutativity in earlier papers (e.g., [^,Q). 

A general form for the Potts model partition function for the strip graphs considered 
here, or more generally, for recursively defined families of graphs comprised of m repeated 
subunits (e.g. the columns of squares of height Ly vertices that are repeated L^ times to 
form an L^ x Ly strip of a regular lattice with some specified boundary conditions), is |pB 



Z{G,q,v)= J2 CgA^gM.vW (2.4) 

i=i 

where the terms Agj, the coefficients cgj, and the total number Nz^g,x depend on G through 
the type of lattice, its width, Ly, and the boundary conditions, but not on the length. 
Following our earlier nomenclature |]73|, we denote a A as leading (= dominant) if it has a 



magnitude greater than or equal to the magnitude of other A's. In the limit n -^ oo the 
leading X in Z determines the free energy per site /. The continuous locus B where / is 
nonanalytic thus occurs where there is a switching of dominant A's in Z and P, respectively. 



and is the solution of the equation of degeneracy in magnitude of these dominant A's. The 
special case of this for the chromatic polynomial was discussed in |7^J75| . 

Let us next define the self-dual strip graphs that we consider here. The first type, 
with boundary conditions that we denote as DBCl (following our nomenclature in |^3|) was 



discussed in Ref. |]38| and illustrated in Fig. 1 of that paper. We shall need a straightforward 
generalization of it to the case of an L^, x Ly lattice with L^ 7^ Ly, and we describe this 
as follows. Let the lattice be oriented with the x and y directions being horizontal and 
vertical, respectively. Let all of the vertices on the upper and right-hand sides, including 
the corner vertices, connect along directions outward from the lattice to a common vertex 
adjoined to this lattice (so that the upper right corner connects to this adjoined vertex via 
edges in both the x and the y directions), while all of the sites on the lower and left-hand 
edges, excluding the previously mentioned corner vertices, have free boundary conditions. It 
is easily checked that this graph is self-dual. Note that it has free longitudinal (horizontal) 
boundary conditions. The number of vertices is n = L^Ly + 1. Graphs of this type were 
used for calculations of Fisher zeros in [0^ , ^] . 

The second type of self-dual strip graph, used in [^ where it was labelled DBC2 [ [75| , 
can be described as follows. Let the L^ x Ly lattice strip have periodic boundary longitudinal 
(=horizontal) boundary conditions and connect all of the vertices on the upper side of the 
strip to a single external vertex, while all of the vertices on the lower side of the strip have 
free boundary conditions |^6|. An illustration of this type of graph is given in Fig. |l|. This 
has also recently been used for calculation of Fisher zeros in [Q. In ||7^ we gave exact results 
for structural properties of Potts model partition functions and chromatic polynomials for 
strips of this type, of arbitrarily great length and width, and presented exact calculations of 
chromatic polynomials and resultant singular loci B for f = — 1 in the q plane for widths up 
to L„ = 4. 
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FIG. 1. Example of an Ly x L^ strip of the square lattice with DBC2 boundary conditions, for the case 
Ly = 3,L,= 4. 



We comment on an interesting feature of the complex-temperature phase diagrams for 
self-dual infinite-length, finite-width strips. In our earlier work yielding exact determinations 
of these complex-temperature phase diagrams for non-self-dual strips [p^, ^ , ^ ,p^-p8| , it was 
found that for some cases, e.g. strips with periodic longitudinal boundary conditions, B 



7 



passes through the origin of the u = a"^ plane. For free longitudinal boundary conditions, 
this does not happen. For the self-dual strip graphs considered here, we can easily prove 
that B does not pass through u = 0, i.e v = ( = oo, since by the inversion symmetry under 
( — > 1/C, this would imply that it passes through C = 0, but this is the infinite-temperature 
{P = 0) point, where the free energy is analytic, so no singular phase boundary can pass 
through this point. 



III. Ly = l STRIP WITH DBCl 

In this section we present the Potts model partition function Z{Ly x Lx, DBCl, q, v) for 
the strips of the square lattice of width Ly and arbitrarily great length L^ = m + 1 containing 
m edges in each horizontal row of the strip, with duality-preserving boundary conditions of 
type 1. We label such a strip graph with DBCl boundary conditions as 5*, Ly or just SLy 
for short and {S, Ly)m to indicate the length. The number of vertices is n = L^Ly + 1. One 
convenient way to express the results is in terms of a generating function, 

oo 

T{S, Ly, q, v,z) = '^ Z{S, Ly, m, q, v)z^ . (3.1) 

m,=0 

As indicated, the coefficients in the Taylor series expansion of this generating function in 
the auxiliary variable z are the partition functions for the strip of length m. We have 
calculated this generating function using transfer matrix methods and iterative application 
of the deletion-contraction theorem for the corresponding Tutte polynomial. We find 

r(S.L,.,....) = ^[fi^^ (3.2) 

V[b,Ly,q,V,z) 

where the numerator Afi^S, Ly, q, v, z) and the denominator V{S, Ly, q, v, z) are polynomials 
in z, q, and v that depend on Ly but not L^. The degree of the denominator in z, i.e., the 



number of A's in the form ( p.4|) , is ||T7 



2 /2/ + l\ 
iVz,5,L„A = deg,{V{S, Ly, q, v, z)) = j-^ ( L ) ' ^^'^^ 

We first treat the minimum-width case, Ly = 1, which has the appeal that the analytic 
results are simple but already exhibit a rich variety of behavior for the locus B. This family 
of graphs can be represented as an open wheel formed by ?7i -|- 1 vertices along the rim, each 
except the rightmost one connected by a spoke (edge) with a vertex forming the axle of 
the open wheel, and with the rightmost vertex on the rim connected by a double edge to 
this central vertex. We calculate for the denominator of the generating function (with the 
abbreviation 5*1 for S,Ly = 1) 

V{S1, q,v,z) = l-{3v + q + v^)z + v{v + l){v + q)z^ 



where 



with 



11(1 - A^i,,^) 



^3,1,(1, 



(1,2) 



Tsi ± x/R 



SI 



Tsi = 3v + q + v'^ 



(3.4) 



(3.5) 



(3.6) 



Rg^ = 5v^ + 2vq + 2v^ + q^ - 2v^q + w^ 

= 5gC' + Sg^/^C + 2g^/'C' + g' - 2g\' + g'C' 
and for the numerator of the generating function 

XiS, 1, q, V, z) = Asifi + Asi,iz 
with 

Asi^o = q{2v + q + v^) 



(3.7) 



(3.^ 



(3.9) 



^51,1 = -qv{v + l){v + q) 



(3.10) 



Ref. I?^ presented a formula to obtain the chromatic polynomial for a recursive family of 
graphs in the form of sums of powers of A^-'s starting from the generating function, and the 



generalization of this to the full Potts model partition function was given in |Q . Using this, 
we have 



Z{Slm,q,v) 



(^51,0-^51,1 + ^51, l) /, ^m , (^51,0^51,2 



A51.2 + A 



(A 



51,1 



^51,2) 



(A 



51,2 



A51,l) 



51,1/ / -V \m 

1^51, 2J 



(3.11) 



It is readily verified that this is symmetric under the interchange Xsi,i ^^ Xsi,2- The free 
energy is given by / = In A^i^i and is analytic for all finite temperature. Regarding the locus 
i3 as a submanifold (with possible singularities) in the C^ space defined by the variables 
(g, v) or (g, (), we can obtain the slices of this locus in the complex q plane for fixed v and 
in the complex f or (^ plane for fixed q. 



1. Ly = 1 with DBCl: Bg for fixed v 

For the physical range v G [—1, oo], the locus B in the q plane consists of a single self- 
conjugate arc which has endpoints at 



qe,q: = viv-l±2tyAi+l) (3.12) 

and crosses the real q axis at 

g = -t,(t, + 3). (3.13) 

These points qe, ql are the branch points of the square root \/Rsi- For the Potts antiferro- 
magnet at T = 0, i.e., v = —1, the locus B degenerates to the single point q = 2. As one 
increases the temperature above T = 0, B expands to form the generic arc as given above, 
but then as T — i> oo, i.e., v — > 0~, this arc shrinks again to a point at the origin, q = 0. 
For the Potts ferromagnet, as T decreases from infinity, i.e., v increases from 0, the arc is 
centered in the negative Re{q) half-plane and crosses the negative real q axis at the value 
given in eq. ( p.l3|) . 



2. Ly = 1 with DBCl: ^^ for fixed q 

We discuss here the locus B in the complex ( plane for fixed q. We shall concentrate on the 
range of real g > 1. The locus B is defined by the equality of magnitudes |A5,i,i| = |As,i,2|- 
This equality can arise in two, in general separate, ways. First, (for real q), on the real 
axis of the ( plane, since Tsi is real, if Rsi < so that the square root in ( |3.5|) is pure 
imaginary, it follows that jA^^i^il = |A5_i^2|- Second, at the four points where Rsi = 0, clearly 
As,i,i = A5_i^2- At certain special values of q some of these six points can coincide. 

Proceeding to analyze B, we first observe that at g = 1, this locus is comprised of a closed 
oval curve that surrounds the point ( = —1 and crosses the real ( axis at C = (l/2)(— 3±-\/5 ), 
i.e., at C ~ —2.618 and —0.3820. The crossing point on the left is at C = —B5, where 

Br = Acos^(-) (3.14) 

is the Tutte-Beraha number. Recall here that B = Bqn, if one were to use the opposite order 
of limits in ( |2.3| ), then Z{G,q = l,v) = {v + 1)" and all zeros collapse to the single point 
V = C = —1. This illustrates the noncommutativity discussed in the introduction. As q 
increases above q = 1, two complex-conjugate arcs sprout out from the points C = e='=^'^*/^, 
so that B is comprised of the union of these arcs and the closed oval curve. The endpoints 
of the arcs occur at the four zeros of Rsi, i-e., branch point zeros of \/Rsi- The right-hand 
endpoints of the arc are located at the complex conjugate pair of points 

1 , ./ . 3\l/2n 



at the angles ±.6ae given by 
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v/^~I--±^(v/^~I+-j J (3.15) 



Oae = arctan 
The left-hand endpoints occur at 

For q < Qac, where 



3 + 4v/^"^T 



1/2 



2v/^"^T-l 



/ 1 / / 3\i/2 



(3.16) 



Qac = (5/4)2 _ ^^5g25 



(3.17) 



(3.18) 



(se is complex and (~J^ = (*^, while for q > (5/4)^, ^^e is real. For q = (5/4)^, we have 
(se = —1- Thus, as q increases through the value q = qac, the left-hand endpoints of the 
complex-conjugate arcs come together and pinch the negative real axis at C = ~1- For 
q > qac, this part of the locus B forms a line segment on the negative real axis centered at 
C = —1, whose right-hand end is (se and left-hand end its inverse. The above-mentioned 
oval curve crosses the negative real ( axis at the two points where Tsi = 0, i.e., 



Ct,C 



-1 



1 



2v/g 



3 ± v^g - 4g 



(3.19) 



For g = 1, Ci = (l/2)(-3 + VS ) ~ -0.3820, C^ = -(V2)(3 + VS ) = -^5, as discussed 
above. The locus B and corresponding complex-temperature phase diagram is plotted for a 
typical value in the interval 1 < g < qac, namely, q = 5/4 = 1.25, in Fig. |[ Here 9ae = 7r/2. 
For comparison, complex-temperature Fisher zeros are shown for a long finite strip. As is 
evident, the density decreases strongly as one approaches the intersection points (multiple 
points in the terminology of algebraic geometry) where the arcs and the oval curves cross 
each other. This is the same behavior that we found in many previous studies of complex- 
temperature zeros for spin models, e.g., [^,^,^,|^. The complex-temperature extension 
of the paramagnetic (PM) phase occupies the full ( plane except for an O region enclosed 
by the oval curve (and the singular set of measure zero comprised by B itself). 
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FIG. 2. Singular locus B in the ( plane for the free energy of the q — 5/4 Potts model on the La, — > cjo 
limit of the Ly ~ 1 strip with DBCl boundary conditions. For comparison, zeros of Z for L^ ~ 21 are 
shown. 

For an interval q > Qac, there are two O phases, namely the regions surrounded by the 
oval curve and separated by the arc of the |C| = 1 circle that passes through ^ = — 1. The 
rest of the ( plane is occupied by the complex-temperature extension of the paramagnetic 
(PM) phase. This type of situation is illustrated in Fig. || for g = (5/4)^ and Fig. ^ for 
g = 1.7. 
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FIG. 3. Same as Fig. |for g = (5/4)2 ^ 15625. 
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FIG. 4. Same as Fig. |for g = 1.7. 

As q increases further, the two O phases that were contiguous now contract and pull 
away from each other. This is illustrated in Fig. ^ for g = 1.85. Eventually, these two O 
phases pull completely away so that they are no longer contiguous; they are then centered 
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around the endpoints of the hne segment. An example is shown in Fig. ^ for q = 1.9. Here 
the crossing of the right-hand O phase, i.e. the crossing nearest to the origin, occurs at Ct 
and its inverse is the leftmost crossing of the other O phase. 
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FIG. 5. Same as Fig. for g = 1.85. 
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FIG. 6. Same as Fig. |for g = 1.9. 
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As q increases to g = 2, the two O phases contract to points and disappear. One can 
see this analytically since at g = 2 the crossing d = — 1/v^ coincides with the line segment 
endpoint Cse = — 1/"\/2 and similarly for their inverses. In the interval 2 < q < 9/4 = 2.25 
the points (t and (f^ are located in the interior of the line segment and do not play a 
special role. As q increases above q = 9/4, there ceases to be any real-g solution of the 
condition T51 = 0. The locus B and complex-temperature phase diagram are shown for 
g = 3 in Fig. |^ A topological feature of B for this region of q, namely the presence of a 
complex-temperature endpoint of a line segment on the left, is reminiscent of the suggestive 
possibility of prongs (or perhaps cusps) on B for the square-lattice Potts model inferred from 
the combination of calculations of Fisher zeros and the correlation of the positions of these 
zeros with locations of complex-temperature singularities that were reliably determined from 
analyses of low-temperature series expansions in |Q (see also [^^ ). The angle 9ae of the 
upper arc endpoint decreases as q increases, i.e., this endpoint moves toward the point ( = 1 
on the real axis. In Table | we list some explicit values of this angle 6ae as a function of q 
for (the La; — >■ 00 limit of) this strip. 



1.5 

1 

0.5 

Im(zeta) 

-0.5 

-1 

-1.5 



T — ' — I — ' — I — ' — I — ' — r 



"-•—*■-•-< 




J , I , I , L 



_J_ 



2 -1.5 -1 -0.5 0.5 
Re(zeta) 



FIG. 7. Same as Fig. | for q = 3. 

TABLE I. Values of arc endpoint angles 9ae for various infinite-length finite-width square lattice strips 
with duality-preserving boundary conditions. The 9ae values are the same for DBCl and DBC2. 



Ly 


q=l 


q = 2 


q = 3 


(7 = 4 


q=W 


1 


120° 


69.3° 


58.1° 


52.0° 


37.8° 


2 


88.8° 


47.7° 


38.8° 


34.1° 


23.7° 
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As q gets large, the right-hand arc endpoints move down toward the point C = 1 ^^^d 
the circular arc finally pinches this point in the limit as g — i> oo. We calculate the following 
expansions for the position of the upper arc endpoint for large q: 

Cae = 1 + iq-''' - \q-"' + l^g-^/^ - \q-' - ^q-'^' + ^^iq-'" - ^^"^ + 0{q-''') 

(3.20) 

(the lower one being the complex conjugate) and, for the right and left endpoints of the line 
segment. On the left, as q gets large the line segment contracts toward C, = —1 and finally 
degenerates to a point at g = oo. We calculate the following expansion for the positions of 
the endpoints of this line segment for large q: 



2^ T^g +-q T^. 



Ce, C/ = -1 ± q-'" - ^q-'" T i,q-"' + -q'^ T 4f^-"' 



^^^<i~"" + Y.<i~" + 0{q~'")- (3-21) 

Thus, in the limit as g — i> oo, i3 becomes the unit circle \C,\ = 1. 

Going the other way, let us start at g = co, where B is the circle |C| = 1. As g decreases 
from infinity, two changes occur immediately in B: (i) the circle breaks open on the right 
side, forming two arcs with endpoints at the angles given in ( ^.16 ) that recede away from the 



real axis, and (ii) a real line segment sprouts out from the point C = — 1- Feature (i) refiects 
the quasi-one-dimensional nature of the L^ —>■ oo limit of this family of strip graphs, since for 
finite g, the free energy of the Potts ferromagnet is analytic for all finite temperatures. This 
feature does not hold in the thermodynamic limit on the square lattice. As was discussed in 
3^,34] for the Ising model and in |^| for the general Potts model, the region around a = 1, 



i.e., V = 0, which is the paramagnetic phase, is not analytically connected to the broken- 
symmetry, ferromagnetic phase; hence the part of the phase boundary B that separates the 
complex-temperature extensions of the PM phase and FM phases from each other must 
remain intact as g decreases from infinity. However, the deviation on the left serves as a 
prototype of the sort of deviations that are suggested by finite-lattice calculations of Fisher 
zeros for sections of the square lattice P^J39|J4^J60|] and gives some insight, as an exactly 



calculable example, of how these deviations arise. As g decreases to sufficiently small values, 
the line segment changes to a region around ( = —1. As we shall show below, the nature of 
the deviation on the left can be more complicated for wider strips. Thus, from the point of 
view of increasing g, B only becomes the unit circle \(\ = 1 in the limit g —>■ oo. 

We next point out an another important feature of these results. For the g = 3 case, B 
crosses the real axis at C = — l/v3, its inverse — v3, and at C = ~1- Transforming back 
to the V plane, these points correspond to f = — 1, f = —3, and v = —VS, respectively. 
The crossing at f = —1, i.e., a = 0, connotes a zero-temperature critical point of the g = 3 
Potts antiferromagnet on this infinite- length Ly = 1 strip p7|j88[] . This is very interesting. 



since this model also has a zero-temperature critical point on the (infinite) square lattice. 
Thus, an infinite-length strip with width Ly = 1 already exhibits a feature of the Potts 
antiferromagnet on the full square lattice. As will be seen below, this is also true of the 
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other widths Ly = 2,3 for which we have obtained exact solutions for the Potts model free 
energy. The crossing at v = —3, i.e., a = —2, is a complex-temperature singular point that 
is the dual image of a = and, by duality, the singularity in the free energy is the same as 
at the physical zero-temperature critical point, in accordance with the general discussion in 



| 45| relating physical and complex-temperature singularities by duality. It is instructive to 
view the locus B in the complex ^ or w plane for fixed g as a slice of the singular subset in 
the full C^ space defined by the pair of variables or (g, v). Thus, the crossing of i3 at f = — 1 
for g = 3 is the point {q,v) = (3, —1) and corresponds to the crossing of the slice of B at 
g = 3 in the g plane at fixed t> = — 1. This is precisely the Qc that we found previously in 
our study of chromatic polynomials and their asymptotic limits and loci B for this family of 



graphs in [77 1 . It will be recalled that we found that this gc = 3 value was universal for all of 



the widths 1 < Lj^ < 4 for which we calculated exact solutions for the chromatic polynomial 
and resultant B. As we shall show below, this corresponds to the feature that for each of 
the widths of strips that we study here, B passes through v = —1, i.e., C = — 1/"\/3 (and, by 
duality, its inverse, ( = — -\/3) for g = 3. 

For the Ising case q = 2, B crosses the v axis at f = — 1 and the dual image v = —2. The 
crossing at f = — 1 connotes a zero-temperature critical point for the Ising antiferromagnet 
on the Lx —>■ oo limit of this graph. 

We have also calculated Fisher zeros for the strips with DBCl and Ly = 2,3. A typical 
example is Ly = 3, g = 4, shown in Fig. |^. For lack of space we do not show the others here, 
but they are available upon request. 
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FIG. 8. Complex-temperature Fisher zeros in the C plane for the q — 4: Potts model on the Ly — 3 strip 
with DBCl boundary conditions with Lx = 21 {n ~ 64). 
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IV. Ly = 1, DBC2 



To elucidate the dependence of B on the self-dual boundary conditions, we consider the 
Ly = 1 strips with DBC2. We shall denote this family generically as L, and L, 1 or LI to 
specify the width. The number of A's in the form (p.4| ), is, from [[77| , 



Nz,L,Ly,X 



2Ly + l 

Ly + l 



In particular, for Ly = 1, this gives A^z,l,i,a = 3. We find 

The corresponding coefficients in ( p.4|) are 

Cli,i = f^^'^ = q{q - 2) 



(4.1) 



(4.2) 



(4.3) 



where 77 



cli, 



M) 



K 



(1) 



q for j = 2, 3 



(4.4) 



f/2d(^)-r2,(^) 



VqU: 



2d-l 



if) 



d-1 
j=0 



.Y2rf-i-j\ .-. 



(4.5) 



where T„(a:) and Un{x) are the Chebyshev polynomials of the first and second kinds. Struc- 
tural properties of Z{G,q,v) for these strips have interesting connections with Temperley- 



Lieb algebras and Bratteli diagrams, which were pointed out in [77 



In Figs. p|-p!3| we show the locus B and associated complex-temperature phase diagram 
in the ( plane for the values g = 2, 3, 4, 5, and 100. (Note that the locus B shown in Fig. 
^ is Bnq-) The arc endpoints of the portion of B lying on the circle |C| = 1 are the same 



as for the [L^ —>■ oo limit of the) Ly = 1 strip with DBCl, i.e., Cae^Ce given in eq. ( |3.15| ). 
The reason for this property is that in this area of the ( plane the locus B is determined by 
the equality in magnitude of two terms, Xl,i,2 and \l,i,3, which are common to the partition 
functions for DBCl and DBC2. In our calculations for Ly = 2,3 we have found the same 
property to hold, so that for a given width Ly and value of q, for the strips considered here, 
the locations of the right-hand arc endpoints are independent of whether one uses DBCl or 
DBC2 boundary conditions. However, there is an interval in q for which A^^i^i is dominant 
in the vicinity of C = —1, and this leads to at least one complex-temperature O phase (in 
the nomenclature of |Q). Figs. W and |TT| illustrate this for the cases g = 3 and g = 4. 
One observes complex-conjugate triple points on B for g = 3. An exactly solved case of a 



1^ 



complex-temperature phase boundary exhibiting such a triple point was given in |42|, where 
it was shown how this triple point results from three curves on B coming together such that 
as one travels along a given curve, beyond the intersection point the A's that were leading 
and degenerate in magnitude on this curve are no longer leading, so that their degeneracy 
is not relevant to B. In our exact calculations of B in the q plane for n — ;> oo limits of 
chromatic polynomials we have found a number of such triple points (e.g., [79|- |84|, flT 



Considering B as the union of the various curves and line segments that comprise it, this is a 
multiple point (intersection point) on B since it lies on multiple branches of B. In a different 
nomenclature in which one considers each of the algebraic curves individually, including the 
portions where the pairs of degenerate-magnitude A's are not dominant so that these portions 
are not on B, then such triple points are not multiple points on each individual algebraic 
curve, since these individual curves pass through the triple point as shown in Fig. 2 in 
(see also 
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5,86 



In general, as we did for the DBCl strips, we find that for any finite q no matter how 
large, B deviates from the circle \(\ = 1. As discussed above, the gap that opens in the circle 
in the vicinity of C = 1 is a property that is special to the quasi-one-dimensional nature 
of these infinite-length, finite width strips. However, just as our previous exact results 
showed that certain complex-temperature properties of quasi-one-dimensional spin models 
were similar to those of the same models in 2D 
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mwm 



so also the deviations in 
the region around ( = —1 are indicative of what can happen in 2D in this case. Note that 
for q > 1, the point ( = —1, i.e., v = —^/q, is a complex-temperature, rather than physical, 
point. We can now use our results to address the issue of the radius of convergence of the 
1/q expansion as regards the form of B. In general, if an expansion of some quantity in 
a variable e has a finite radius of convergence ec, then, roughly speaking, the behavior for 
|e| < ec should be qualitatively the same as for e = 0. Our results suggest that, at least 
for the infinite-length finite-width strips, the 1/q expansion has zero radius of convergence 
insofar as properties of the locus B are concerned. This follows since B in the ( plane 
differs qualitatively for any nonzero value of 1/q from its form at 1/q = 0. We know 
that the deviation near ^ = 1, i.e. the physical PM-FM transition temperature given by 
Vc = ^/q, that occurs for these quasi-one-dimensional strips will be absent in 2D. However, 
the deviation in the vicinity of the complex-temperature point ( = —1 should be a more 
general feature, not limited to the quasi-one-dimensional nature of the strips considered 
here. This inference follows from (i) our previous experience comparing exactly determined 
complex-temperature features of Ising model phase diagrams for infinite-length strips and 
in 2D, (ii) the observed scatter of Fisher zeros in 2D ||29|j38|j43[| in this complex-temperature 
region, and (iii) reliable determinations of locations of complex-temperature singularities 



via analyses of low-temperature series ||43| , ^ , ^ , ^ and the correlation of these with points 
on B ||33|j37|j43|j47|j48|l . Hence for the locus B, our present exact results suggest that the 
1/q expansion has zero radius of convergence. We emphasize that this does not reduce the 
value of these large-g expansions, since the point where the deviation occurs is generically 
a complex-temperature point, and this type of deviation does not occur near the physical 
PM-FM phase transition point. Indeed, large-g expansions yield excellent agreement with 
Monte Carlo measurements of thermodynamic quantities in the Potts model p2||5^ and also 
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for the ground state degeneracy W{q) in the T = Potts antiferro magnet |73 



eel- [72 



Furthermore, it is easily checked that the dominant term A that determines the free energy 
on the strips that we consider does have a well-defined large-g expansion (in the variable 
1/y/q). For example, for Ly = 1 with DBCl or DBC2, removing the leading factor of q, we 
have the expansion for q^^Xs,i,i for large q: 



Q ^^sxi 






C^ 



-Q 



-1/2 



C^(C^ + C-i)( C^-C-i) -1 



iC 



C^(C^ + i)(C^ + C-i)(C^-C-i) -3/2 



0{q 



(4.6) 



Note the poles aX C, = ±1 in this expansion. Parenthetically, we note that our findings here 
concerning the large-g expansion are not related to our earlier studies of families of graphs 
for which W{q) has no large-g expansion [ ]75| , pU| , pD| , pT| , ^ (see also ||5T|) since in those cases, 
the breakdown of the 1/g expansion is equivalent to the property that the singular locus 
B is noncompact in the q plane, passing through the origin of the 1/q plane. This sort of 
breakdown does not occur for the present families of graphs, as is clear from the fact that B 
is compact in the q plane, shown in [^ and above. 
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FIG. 9. Singular locus B in the C plane for the free energy of the q — 2 Potts model on the L^; — > oo limit 
of the Ly — 1 strip DBC2 boundary conditions. For comparison, zeros of Z for Lx ~ 20 are shown. 
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FIG. 10. Same as Fig. | for g = 3. 
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FIG. 11. Same as Fig. | for <? = 4. 
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FIG. 12. Same as Fig. | for g = 5. 
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FIG. 13. Same as Fig. | for g = 100. 
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V. WIDER STRIPS 



We have also calculated Z{G, q, v) for arbitrary q and v for wider strips with Ly = 2 and 



Ly = 3, for which our general formulas ( p.3| ) and ( [4.1| ) from |^ yield for the number of A's 
the results Nz,dbci,2,x = 5, Nz,dbc2,2,x = 10, Nz,dbci,3,x = 14, and Nz,dbc2,3,x = 35. The 
analytic expressions for the A's are too lengthy to list here. In several figures below we show 
plots of Fisher zeros in the ( plane for various values of q. 
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FIG. 14. Complex-temperature Fisher zeros in the C plane for the q — 2 Potts model on the Ly = 2 strip 
with DBC2 boundary conditions with L^ = 20 {n — 40). 
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FIG. 15. Same as Fig. plfor g = 3. 
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FIG. 16. Same as Fig. ||for g = 4. 
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FIG. 17. Same as Fig. || for g = 100. 
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FIG. 18. Complex-temperature Fisher zeros in the C plane for the q = 2 Potts model on the Ly = 3 strip 
with DBC2 boundary conditions with L^. = 20 (n = 60). 
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FIG. 19. Same as Fig. |9|for g = 3. 
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FIG. 21. Same as Fig. [| for q = 100. 

It is interesting to note that the simplification of the complex-temperature phase diagram 
for the Potts model to the circle |C| = 1 for q ^ oo proved in ^§\ and studied further here 
with exact results is somewhat similar to the simplification of the complex-temperature phase 
diagram for the 2D spin-s Ising model in the limit s ^ cxd to the circle \us\ = 1 in the plane 
of the Boltzmann variable Us = e~^l^ p^j93|J^ . In both cases, the approach to the limit is 



singular in the sense that there are deviations from the asymptotic locus. 



VI. CONCLUSIONS 

In summary, we have presented exact calculations of the Potts model partition function 
Z[G^ g, v) on self-dual strip graphs G of the square lattice with fixed width Ly and arbitrarily 
great length L^ with two types of boundary conditions. In the infinite-length limit we have 
studied the resultant complex-temperature phase diagram. In particular, we have considered 
the widths Ly = 1, 2, 3. We have used these results to study the approach to the large-g 
limit of B, where this locus is the unit circle \(\ = 1. We find that, for a given Ly and set of 
self-dual boundary conditions, a portion of B lies on this unit circle, while for any finite q, a 
portion deviates from this circle. For the strips considered here we find that the right-hand 
arc endpoints on the circle are independent of whether the longitudinal boundary conditions 
are periodic or free and to curve around and pinch the real axis at ^ = 1 as g — >■ cxd. For a 
fixed value of q, as Ly increases, the right-hand arc endpoints move closer to C = 1- On the 
left, the nature of the complex-temperature phase diagram was found to depend in detail on 
both the type of boundary conditions and the width of the strip. As g ^ oo, the deviations 
typically include real line segments as well as possible O phases. One feature was found 
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for each width considered, namely that for the DBC2 strips, for q = 3, B crosses the ( at 
( = —l/\/3 and at — v^- We showed that this is equivalent to the fact that for f = —1, each 
of the infinite-length strips that we studied in |j7^ with I < Ly < 4 and DBC2 had gc = 3 
as for the infinite square lattice. As discussed, the gap in the locus B that opens around 
C = 1 as 1/g increases above zero is a consequence of the quasi-one-dimensional nature of 
the strips. However, the behavior in the Re{() < region near ( = —1 can give some insight 
as to how B could behave for large q on the infinite square lattice. 

Acknowledgment: This research was partially supported by the NSF grant PHY-97- 
22101. 

VII. APPENDIX 

A. General 

The most compact way to express a Potts model partition function Z{G, q, v) is often in 
terms of the corresponding Tutte polynomial T{G,x,y). For the reader's convenience, we 
recall the definition of the Tutte polynomial and some basic formulas relating these functions 



here (e.g., ^3[). For an arbitrary graph G the Tutte polynomial of G, T(G, x, y), is given by 

i-0 

r(G,x,y)= 5^(x-i)^(«')-'=(^)(y-ir(^') (7.1) 

G'CG 

where the spanning subgraph G' was defined in the introduction, and we recall that k{G'), 
e{G'), and n{G') = n{G) denote the number of components, edges, and vertices of G', where 

c{G') = e{G') + k{G') - n{G') (7.2) 

is the number of independent circuits in G'. As stated in the text, k{G) = 1 for the graphs 
of interest here. Now let 

X = 1 -\ — , y = a = V + 1 (7.3) 

V 

so that 

q = {x-l){y-l) . (7.4) 

Then 

Z{G, q, v) = {x- l)^(^) {y - l)"(«)r(G, X, y) . (7.5) 

For a planar graph G the Tutte polynomial satisfies the duality relation 

T{G,x,y) = T{G*,y,x) (7.6) 



where G* is the (planar) dual to G. As discussed in ||S^ , the Tutte polynomial for recursively 
defined graphs comprised of m repetitions of some subgraph has the form 
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Nx 






(7.7) 



One special case of the Tutte polynomial of particular interest is the chromatic polynomial 
P{G, q). This is obtained by setting y = 0, i.e., v = —1, so that x = 1 — g; the correspondence 
is P{G, q) = (-g)'=(^)(-l)"T(G', 1 - g, 0). 

B. Strips with DBCl 

The generating function representation for the Tutte polynomial for the strip Sm of the 
square lattice with length L^ = m + 1 vertices, i.e., m edges in each horizontal row, and of 
width Ly, with duality preserving boundary conditions of type 1, is 



^TWmi Ly, X, y, Z) — y ^ T[Sm, Ly, X, %)) 



z^ . 



m=0 



We have 



VT{S,Ly,x,y\z) 



MT{S,Ly,x,y,z) 
VTiS,Ly,x,y,z) 



For L,, = 1 we find 



(7i 



(7.9) 



Mt{S, 1, X, y, z) = (x + y) - xyz 



(7.10) 



Vt^S, l,x,y,z) = 1 — {1 + X + y)z + xyz^ 



with 



At,5,1,(1,2) - - 



n(l - At,5,1,,^) 



l + x + y±(l + 2{x + y) + {x- yf\ 



1/2 



(7.11) 



(7.12) 



T{S„i,x,y) 



^T,S,O^T,S,l + ^T, 



5,1 



At,5,1 — At,5,2 



(At.s'.i)™' + 



At,S,0^T,S,2 



At,s,i 



^T,S,2 — ^T,S,1 



{\T,s,2r ■ (7.13) 



For Ly = 2 we find 



AfriS, 2, X, y, z) = As2fl + As2,iz + ^52,2^' + ^52,3^^ + ^52,4^' 



where 



(7.14) 



2 I „.2 



As2,o = x + y + xy + x +y 



(7.15) 
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As2,l = — 



x + y + 2{x + y ) + 3xy + 5xy{x + y) + xy{x +y) + x +y + {xy) (7.16) 



^52,2 = xy 3{x + y)+ 4(a;^ + y^) + 6xy + 3xy{x + y) + x^^ + y^ + {xyY (7.17) 



As2,3 = -ixy) 2{x + y)+x +y +3xy + xy{x + y) 



{7.U 



^52,4 = {xyY 



(7.19) 



Vt{S, 2, X,y,z) = 1 + bs2,lZ + bs2,2Z^ + bs2,3Z^ + bs2,4Z^ + bs2,5Z^ 



where 



2 I „,2 



^■52,1 = ~ 3{1 + X + y) + xy + X +y 



(7.20) 



(7.21) 



&52,2 = l + 3{x + y) + 3{x^ + y') + 8xy + {x'^ + y^) + 5xy{x + y) + xy{x^ + y^) + {xyY 

(7.22) 



^52,3 = -xy 3 + 5{x + y)+ 4(x^ + y^) + Qxy + (x'^ + y^) + 3xy{x + y) + {xyy (7.23) 



bs2A = (xyfil + x){l+y)il + x + y) 



(7.24) 



^52,5 



■{xyY ■ 



(7.25) 



For Ly = 3 our general results in [|7^] yield the result that there are 14 terms, and we 
find that the At,5,3,j's are roots of an algebraic equation of degree 14. This is too lengthy to 
record here, but is available upon request. 

C. Strips with DBC2 



We have 



^T.L,Ly.X 

T{Lm,Ly,X,y) = ^ CT,L,Ly,j{>^T,L,L, 



y^3' 



(7.26) 



where NT,L,Ly,x = Nz.L.Ly,x, and our general formula for this number, from |^ was given in 
XTH 
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l.Ly = l 



We have Nt l i a = 3 and 



A 



T,L,1,1 



(7.27) 



AT,L,i,i = At,5,ij-i for J =2, 3 (7.28) 

where \t,s,i,i and \t,s,i,2 were given in eq. (|7.12|) . The corresponding coefficients are 

(7.29) 



Ct,l,i,i = q ^ft:^^^ =q-2 = xy-x-y-l 



CT,L- = q-^n^^) = l for j = 2,3 
where k^'^' was defined in eq. ( [4.5|) in the text. 



/i/. i/y — .z 



We have A^t,l,2,a = 10. The term At,l,2,i is 

Ar,L,2,i = 1 • 
The Xt,l,2,j for 2 < j < 5 are solutions to the equation 



e 



2{x + y) + 3 ^^ + x^ + y^ + 3{x + y)+ Axy + 1 



e 



-xy 2{x + y) + 3 ^ + (xy)"^ = . 
The Xt,l,2,j for 6 < j < 10 are solutions to the equation 

e' + ^52,1^' + ^52,2^' + ^52,3^' + ^52,4^ + &52,5 = . 

The corresponding coefficients are 

Ct,l,2,i = g^^K^^^ = (q - l){q - 3) = {xy - X - y){xy - X - y - 2) 

CT,L,2,j = q'^K^'^^ =q-2 = xy-x-y-l for 2<j<5 



Cr,L,2j = g^^K^^^ = 1 for 6 < j < 10 . 



(7.30) 



(7.31) 



(7.32) 

(7.33) 

(7.34) 
(7.35) 
(7.36) 
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O. Ijy O 

Here our general formulas in |]7^ yield the results that there are 35 terms in all, comprised 
of (i) one term with coefficient g~^K*^^\ namely, At,l,3,i = 1, (ii) six terms At,l,3j, 2 < j < 7, 
with coefficient q~^tv~'^\ (iii) 14 terms \t,l,3,j for 8 < j < 21 with coefficient q~^K^'^\ and (iv) 
14 terms Xt,l,3,j with 22 < j < 35 with coefficient q~^K,^^'. The terms in (ii) are roots of the 
sixth-degree equation 

-(x^ + 9yx2 + 5x2 _j_ g^2^ _j_ 2Qyx ^Qx + y^ + 5y2 + 6y + 1)^^ 

+xy{3x^ + %x + lOx + 3y2 + lOy + 6)^^ _ x^y^(^^x + 3y + 5)^ + x^y^ = . (7.37) 

The equation of degree 14 for the \t,l,3,j with coefficient q~^K^^^ is the same as the single 
degree-14 equation for the terms in the Ly = 3 strip with DBCl. Both this and the other 
degree-14 equation are too lengthy to list here, but can be provided at request. 

D. Special Values of Tutte Polynomials 

For a given graph G = (y,E), at certain special values of the arguments x and y, the 
Tutte polynomial T{G,x,y) yields quantities of basic graph-theoretic interest [[7|- ||TD|. We 



recall some definitions: a spanning subgraph was defined at the beginning of the paper; a 
tree is a connected graph with no cycles; a forest is a graph containing one or more trees; 
and a spanning tree is a spanning subgraph that is a tree. We recall that the graphs G that 
we consider are connected. Then the number of spanning trees of G, Nst{G), is 

Nst{G)=T{G, 1,1), (7.38) 

the number of spanning forests of G, Nsf{G), is 

Nsf{G)=T{G, 2,1), (7.39) 

the number of connected spanning subgraphs of G, Ncssg{G), is 

Ncssg{G)=T{G, 1,2), (7.40) 

and the number of spanning subgraphs of G, Nssg{G), is 

Nssg{G)=T{G, 2,2). (7.41) 



Since the graphs that we consider are self-dual, and using the symmetry relation (|7.6| ), we 
have 

NsFiG) = Ncssg{G*) = Ncssg{G) . (7.42) 

From our calculations of Tutte polynomials, we find the following results. 
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E. Ly = 1, DBCl 



(7.44) 



NsF{Slm) = NcssG{Slm) = (^- + vi\{2 + ^2 T + {^-- vi\{2 - ^2 



NssciSlm) = 2'^^^™^ = 22("^+i) . (7.45) 

Since these quantities grow exponentially, it is natural to define an associated quantity 
that measures this growth [p^ -p6| . In particular, for the number of spanning trees, we define 

Z{G} = lim n"^ In NsT{{G}) ■ (7.46) 

n — ^oo 

For the present Ly = 1, DBCl strips, we thus have z = ln[(3 + -\/5 )/2] ~ 0.9624. A general 
upper bound on the number of spanning trees of a graph G is [^ 

Nst{G)<-(^) , (7.47) 

For the present Ly = 1, DBCl strips, this gives the upper bound 2; < 21n2 ~ 1.386, which 
is seen to be satisfied by our result. 



F. Ly = 1, DBC2 

NsHLl^) = -2 + ('-±^ Y + ('-±^ ) (7.48) 



NsFiLlm) = NcssciLlm) = -2 + {2 + V2r + {2-V2r (7-49) 

NssaiLl^) = 2<'^'-^ = 2^- . (7.50) 

As discussed before |^^, for a given Ly and set of transverse boundary conditions, the 
value of z is the same, independent of whether the longitudinal boundary conditions are free, 
as in DBCl or periodic, as for DBC2. 
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G. Ly = 2, DBC2 

5 10 

Nst{L2^) = 3 - 2 5^(At,l,2,,(1, 1) T + $^(AT,L,2,i(l, 1) T (7.51) 

j=2 i=6 

where At,l,2,j(1, 1), 2 < j < 5 and 6 < j < 10 are the roots of eq. ( [7. 321) and eq. (|7.33| ), 
respectively,' for x = l,y = 1, viz., ^^ - 7^^ + 13^2 _ 7^ + ^ = g, and (^ - 1)(^^ - 11^^ ^ 
25^2- lie +1)=0. 



Nsf{L2„,) = Ncssg{L2J = 3-2 



'3 + v/5\"^ /3-V5 



(l±-)%(i-^)%,3H-v^,. 



10 

+ (3-v^ 



where At,l,2j(2, 1), 6 < j < 10 are the roots of eq. ( [7.33| ) for a; = 2, ?/ = 1, viz., ^^ — 19^^ + 
94^3 _ 162^2 ^ gg^ - 16 = 0. 

NssGiL2„.) = 2^^^''"^ = 2^" . (7.53) 

Hence, in particular, for spanning trees, we find z ~ 1.044 for Ly = 2. It is straightforward to 
use our exact calculations of Tutte polynomials for Ly = 3 with DBCl and DBC2 boundary 
conditions to list similar results for spanning trees, etc. 
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